TABLE XXII. Lattice and experimental determinations of the mass differences of the charmed 
heavy-light mesons (anig = oo). 



mass difference 


lattice (MeV) 


expt. (MeV) 


M(D*) - M(D°) 


38(3) 


146 


M{Di) - [M{D^) + ■iM{D*)]lA 


670(60) 


450 
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TABLE XVIII. Lattice and experimental determinations of the masses of the charmed 
heavy-light mesons (anig = 0.010). 



meson 




lattice (MeV) 


expt. (MeV) 






i700(i00j 




D* 




1750(100) 


2010 


Di 




ZZOU^ lOU J 




TABLE XIX. Lattice 


and experimental determinations 


of the masses of the charmed 


heavy-light mesons (anig = 


0.025). 




meson 




lattice (MeV) 


expt. (MeV) 






1 RQCif 1 1 n"\ 

iDyu(^iiu ) 


ioD4 


D* 




1730(110) 


2010 


Di 




zuyu^ i^u J 


z^zo 


TABLE XX. Lattice an 


d experimental determinations of the 


masses of the charmed heavy-light 


mesons (anig = oo). 






meson 




lattice (MeV) 


expt. (MeV) 






1670(140) 


1864 


D* 




1740(140) 


2010 






2050(170) 


2423 



TABLE XXI. Lattice and experimental determinations of the mass differences of the charmed 
heavy-light mesons (anig = 0.010). 



mass difference 


lattice (MeV) 


expt. (MeV) 


M(D*) - M(D°) 


50(3) 


146 


M{Di) - [M{D^) 3M(i:i*)]/4 


560(30) 


450 
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TABLE XV. Mass Splittings and inverse lattice spacings for the anig = 0.025 configurations. 



Mass difference 


AM (MeV) 


K 


aAM 


(MeV) 


Using tp^astp 










aAMi^Pi - S) 


442 


0.f320 
0.f4fO 


0.257(f0) 
0.274(ff) 


f660(40) 


Using cubic group 










aAM(3Pi - S) 


442 


0.f320 
0.f4fO 


0.256(2f) 
0.302(36) 


f580(f20) 


aAM^Pi - S) 


458 


0.f320 
0.f4fO 


0.249(f7) 
0.284(25) 


f720(f00) 


aAM(P - S) 


457 


0.f320 
0.f4fO 


0.258(26) 
0.293(26) 


f660(ff0) 



Overall estimate of a"^ f660 ± f fO (stat.) ±fOO (sys.) MeV 



TABLE XVL Inverse lattice spacings (in MeV) using other quantities. 



aniq 


"force" 


m,(W) 




mp(S) 


mp(S) 


NRQCD 


0.025 


1935 


2000 


1685 








0.010 


2055 


2140 


1800 






2400 


0.0 


2135 


2230 


1875 


1800 


1660 





TABLE XVIL Mass Splittings and inverse lattice spacings for the quenched (anig = oo) con- 



figurations. 



Mass difference 


AM (MeV) 


K 


aAM 


(MeV) 


aAM(3Pi - S) 


442 


0.1300 
0.1450 


0.23(3) 
0.26(3) 


1800(150) 



Overall estimate of a"! 1800 ± 150 (stat.) ±100 (sys.) MeV 
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TABLE XIV. Mass Splittings and inverse lattice spacings for the anig = 0.010 configurations. 



Mass difference 


AM (MeV) 


K 


aAM 


(MeV) 


Using tp^astp 










aAMi^Pi - S) 


442 


0.1320 
0.1410 


0.238(9) 
0.248(9) 


1810(50) 


Using cubic group 










aAM(3Pi - S) 


442 


0.1320 
0.1410 


0.235(16) 
0.248(20) 


1800(100) 


aAM^Pi - S) 


458 


0.1320 
0.1410 


0.234(16) 
0.250(18) 


1900(90) 


aAM(P - S) 


457 


0.1320 
0.1410 


0.234(6) 
0.247(7) 


1900(30) 



Overall estimate of a"^ 1900 ± 50 (stat.) ±100 (sys.) MeV 
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TABLE XIII. Best range fits to M("^Pi) for the quenched data using the axial vector current. 



i 


i 






error 




conf. Ivl. 




h 








0.1300 


1.650 


0.019 


9.42/13 


0.583 


10 


22 


1 





0.1375 


1.376 


0.025 


6.80/13 


0.815 


10 


22 


1 


1 


0.1450 


1.093 


0.038 


8.27/13 


0.689 


10 


22 


2 





0.1410 


1.237 


0.035 


6.63/13 


0.828 


10 


22 


2 


1 


0.1485 


1.072 


0.024 


13.49/16 


0.489 


7 


22 


2 


2 


0.1520 


0.975 


0.037 


6.43/10 


0.599 


7 


16 


3 





0.1420 


1.317 


0.021 


13.91/16 


0.456 


7 


22 


3 


1 


0.1495 


1.065 


0.029 


14.36/14 


0.278 


7 


20 


3 


2 


0.1530 


0.972 


0.046 


6.59/10 


0.582 


7 


16 


3 


3 


0.1540 


0.967 


0.057 


8.39/10 


0.396 


7 


16 


4 





0.1425 


1.315 


0.023 


11.37/16 


0.657 


7 


22 


4 


1 


0.1500 


1.059 


0.032 


7.83/9 


0.348 


7 


15 


4 


2 


0.1535 


0.986 


0.052 


8.31/10 


0.404 


7 


16 


4 


3 


0.1545 


0.974 


0.062 


10.85/10 


0.210 


7 


16 


4 


4 


0.1550 


0.956 


0.070 


4.35/6 


0.360 


7 


12 
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TABLE XII. Best range fits to M('^Si) for the quenched data using the vector current. 



i 


i 


(Kj + Kj)/2 


aM(^S\) 


error 




conf. Ivl. 




h 








0.1300 


1.425 


0.002 


3.76/6 


0.439 


12 


17 


1 





0.1375 


1.137 


0.002 


7.47/9 


0.381 


14 


22 


1 


1 


0.1450 


0.836 


0.002 


6.02/8 


0.421 


15 


22 


2 





0.1410 


1.003 


0.003 


6.19/9 


0.518 


14 


22 


2 


1 


0.1485 


0.700 


0.003 


4.78/7 


0.444 


11 


17 


2 


2 


0.1520 


0.551 


0.004 


4.41/7 


0.492 


11 


17 


3 





0.1420 


0.971 


0.003 


6.95/9 


0.434 


11 


19 


3 


1 


0.1495 


0.660 


0.004 


2.94/7 


0.709 


11 


17 


3 


2 


0.1530 


0.505 


0.005 


3.85/7 


0.571 


12 


18 


3 


3 


0.1540 


0.460 


0.007 


4.68/8 


0.585 


12 


19 


4 





0.1425 


0.953 


0.004 


3.50/9 


0.835 


11 


19 


4 


1 


0.1500 


0.639 


0.004 


3.43/8 


0.754 


11 


18 


4 


2 


0.1535 


0.483 


0.007 


4.82/8 


0.568 


12 


19 


4 


3 


0.1545 


0.435 


0.009 


4.47/8 


0.614 


12 


19 


4 


4 


0.1550 


0.406 


0.009 


7.96/10 


0.437 


12 


21 
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TABLE XI. Best range fits to M(^So) for the quenched data using the pseudoscalar current. 



i 


i 


(Kj + Kj)/2 


aM(\So) 


error 




conf. Ivl. 




h 








0.1300 


1.411 


0.002 


0.88/5 


0.829 


17 


21 


1 





0.1375 


1.119 


0.002 


1.52/6 


0.822 


17 


22 


1 


1 


0.1450 


0.804 


0.002 


2.06/6 


0.724 


17 


22 


2 





0.1410 


0.978 


0.003 


0.58/6 


0.965 


17 


22 


2 


1 


0.1485 


0.650 


0.002 


5.34/8 


0.501 


15 


22 


2 


2 


0.1520 


0.479 


0.002 


4.30/8 


0.636 


10 


17 


3 





0.1420 


0.937 


0.003 


1.36/6 


0.852 


17 


22 


3 


1 


0.1495 


0.607 


0.002 


7.73/10 


0.460 


12 


21 


3 


2 


0.1530 


0.421 


0.002 


10.57/14 


0.566 


9 


22 


3 


3 


0.1540 


0.364 


0.003 


11.31/14 


0.503 


9 


22 


4 





0.1425 


0.917 


0.004 


2.55/6 


0.637 


17 


22 


4 


1 


0.1500 


0.586 


0.002 


8.23/12 


0.606 


10 


21 


4 


2 


0.1535 


0.393 


0.002 


8.61/14 


0.736 


9 


22 


4 


3 


0.1545 


0.331 


0.002 


8.96/14 


0.707 


9 


22 


4 


4 


0.1550 


0.298 


0.003 


10.68/15 


0.638 


8 


22 
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TABLE X 


. Best ran^ 


le fits to M(iPi) 


for aniq = 


0.025 using the cubic group ^Pi 


state. 




i 


J ( 




aM(iPi) 


error 


xVNdf 


conf. Ivl. 


t^ 




1 


1 


0.1320 


1.715 


0.025 


0.75/3 


0.860 


8 


12 


2 


1 


0.1365 


1.564 


0.027 


1.75/3 


0.627 


8 


12 


2 


2 


0.1410 


1.401 


0.026 


1.99/3 


0.574 


8 


12 


3 


1 


0.1422 


1.371 


0.029 


3.95/3 


0.267 


8 


12 


3 


2 


0.1467 


1.207 


0.029 


4.15/3 


0.246 


8 


12 


3 


3 


0.1525 


0.964 


0.031 


2.24/3 


0.525 


8 


12 


4 


1 


0.1442 


1.296 


0.034 


4.02/3 


0.260 


8 


12 


4 


2 


0.1487 


1.139 


0.036 


3.94/3 


0.268 


8 


12 


4 


3 


0.1545 


0.898 


0.043 


2.22/3 


0.528 


8 


12 


4 


4 


0.1565 


0.804 


0.052 


1.41/3 


0.703 


8 


12 


5 


1 


0.1452 


1.321 


0.021 


12.44/11 


0.331 


4 


16 


5 


2 


0.1497 


1.164 


0.019 


11.61/11 


0.394 


4 


16 


5 


3 


0.1555 


0.855 


0.051 


7.31/7 


0.398 


8 


16 


5 


4 


0.1575 


0.872 


0.037 


12.22/9 


0.201 


6 


16 


5 


5 


0.1585 


0.867 


0.025 


15.33/5 


0.009 


4 


10 


6 


1 


0.1460 


1.163 


0.049 


5.19/7 


0.637 


8 


16 


6 


2 


0.1505 


1.125 


0.024 


9.36/5 


0.096 


4 


10 


6 


3 


0.1562 


0.909 


0.044 


3.48/3 


0.324 


6 


10 


6 


4 


0.1583 


0.851 


0.052 


3.51/3 


0.320 


6 


10 


6 


5 


0.1593 


0.893 


0.121 


3.68/2 


0.159 


7 


10 


6 


6 


0.1600 


0.815 


0.039 


13.27/5 


0.021 


4 


10 
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TABLE IX 


. Best ran^ 


fits to M(3Po) 


for aniq = 


0.025 using 


the cubic group Ai 


state. 




i 


j (k, + Kj)/2 


aM(3Po) 


error 


xVNdf 


conf. Ivl. 


t^ 




1 


1 


0.1320 


1.708 


0.030 


0.28/3 


0.964 


8 


12 


2 


1 


0.1365 


1.550 


0.033 


0.19/3 


0.979 


8 


12 


2 


2 


0.1410 


1.375 


0.034 


0.07/3 


0.996 


8 


12 


3 


1 


0.1422 


1.338 


0.041 


1.65/3 


0.649 


8 


12 


3 


2 


0.1467 


1.154 


0.039 


1.63/3 


0.654 


8 


12 


3 


3 


0.1525 


0.876 


0.040 


1.03/3 


0.794 


8 


12 


4 


1 


0.1442 


1.242 


0.049 


4.98/3 


0.173 


8 


12 


4 


2 


0.1487 


1.057 


0.046 


5.35/3 


0.148 


8 


12 


4 


3 


0.1545 


0.776 


0.047 


3.43/3 


0.330 


8 


12 


4 


4 


0.1565 


0.715 


0.077 


1.83/2 


0.401 


9 


12 


5 


1 


0.1452 


1.186 


0.062 


8.09/3 


0.044 


8 


12 


5 


2 


0.1497 


1.051 


0.218 


9.20/3 


0.027 


8 


12 


5 


3 


0.1555 


1.168 


0.025 


14.69/5 


0.012 


10 


16 


5 


4 


0.1575 


0.935 


0.032 


14.56/6 


0.024 


9 


16 


5 


5 


0.1585 


0.853 


0.032 


16.93/5 


0.005 


4 


10 


6 


1 


0.1460 


1.138 


0.071 


8.63/3 


0.035 


8 


12 


6 


2 


0.1505 


1.005 


0.048 


9.84/3 


0.020 


8 


12 


6 


3 


0.1562 


0.925 


0.029 


18.35/5 


0.003 


4 


10 


6 


4 


0.1583 


0.856 


0.033 


18.97/5 


0.002 


4 


10 


6 


5 


0.1593 


0.824 


0.036 


17.19/5 


0.004 


4 


10 


6 


6 


0.1600 


0.806 


0.045 


12.93/5 


0.024 


4 


10 
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TABLE VIII. Best ran^ 


^e fits to M(3Pi) 


for aniq 


= 0.025 using 


the cubic group 


Tia state. 




i 


i 






error 


xVNdf 


conf. Ivl. 


t^ 




1 


1 


0.1320 


1.714 


0.028 


0.82/3 


0.846 


8 


12 


2 


1 


0.1365 


1.563 


0.031 


1.60/3 


0.660 


8 


12 


2 


2 


0.1410 


1.398 


0.031 


1.69/3 


0.639 


8 


12 


3 


1 


0.1422 


1.370 


0.038 


3.28/3 


0.351 


8 


12 


3 


2 


0.1467 


1.208 


0.039 


3.36/3 


0.339 


8 


12 


3 


3 


0.1525 


0.994 


0.050 


2.51/3 


0.473 


8 


12 


4 


1 


0.1442 


1.303 


0.046 


4.18/3 


0.242 


8 


12 


4 


2 


0.1487 


1.152 


0.054 


3.99/3 


0.263 


8 


12 


4 


3 


0.1545 


1.077 


0.073 


4.36/4 


0.360 


7 


12 


4 


4 


0.1565 


0.956 


0.041 


2.41/2 


0.300 


9 


12 


5 


1 


0.1452 


1.284 


0.067 


4.29/3 


0.232 


8 


12 


5 


2 


0.1497 


1.260 


0.063 


6.21/4 


0.184 


7 


12 


5 


3 


0.1555 


0.985 


0.020 


7.10/5 


0.213 


4 


10 


5 


4 


0.1575 


0.921 


0.023 


8.47/5 


0.132 


4 


10 


5 


5 


0.1585 


0.894 


0.027 


9.67/5 


0.085 


4 


10 


6 


1 


0.1460 


1.433 


0.112 


5.90/4 


0.207 


7 


12 


6 


2 


0.1505 


1.331 


0.113 


5.95/4 


0.203 


7 


12 


6 


3 


0.1562 


0.951 


0.025 


10.12/5 


0.072 


4 


10 


6 


4 


0.1583 


0.889 


0.029 


12.36/5 


0.030 


4 


10 


6 


5 


0.1593 


0.861 


0.035 


14.33/5 


0.014 


4 


10 


6 


6 


0.1600 


0.837 


0.044 


14.27/5 


0.014 


4 


10 
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TABLE VII. Best ran^ 


fits to M(3P2) 


for aniq 


= 0.025 using 


the cubic group Ea 


state. 




i 


i 




aM(3P2) 


error 


xVNdf 


conf. Ivl. 


t^ 




1 


1 


0.1320 


1.710 


0.026 


1.97/3 


0.579 


8 


12 


2 


1 


0.1365 


1.557 


0.029 


3.68/3 


0.298 


8 


12 


2 


2 


0.1410 


1.392 


0.029 


4.09/3 


0.252 


8 


12 


3 


1 


0.1422 


1.388 


0.033 


8.51/4 


0.075 


7 


12 


3 


2 


0.1467 


1.219 


0.034 


7.83/4 


0.098 


7 


12 


3 


3 


0.1525 


1.018 


0.045 


7.31/4 


0.120 


7 


12 


4 


1 


0.1442 


1.363 


0.063 


9.44/3 


0.024 


8 


12 


4 


2 


0.1487 


1.124 


0.084 


2.79/2 


0.248 


9 


12 


4 


3 


0.1545 


1.035 


0.070 


8.04/4 


0.090 


7 


12 


4 


4 


0.1565 


7.799 


0.000 


22.25/2 


0.000 


9 


12 


5 


1 


0.1452 


1.344 


0.092 


2.64/2 


0.267 


9 


12 


5 


2 


0.1497 


1.242 


0.033 


3.89/2 


0.143 


9 


12 


5 


3 


0.1555 


0.938 


0.027 


8.95/5 


0.111 


4 


10 


5 


4 


0.1575 


0.879 


0.034 


11.12/5 


0.049 


4 


10 


5 


5 


0.1585 


0.866 


0.042 


11.74/5 


0.038 


4 


10 


6 


1 


0.1460 


1.085 


0.114 


0.53/2 


0.768 


9 


12 


6 


2 


0.1505 


0.888 


0.077 


0.81/3 


0.847 


8 


12 


6 


3 


0.1562 


0.908 


0.035 


4.35/5 


0.500 


4 


10 


6 


4 


0.1583 


0.849 


0.042 


5.53/5 


0.354 


4 


10 


6 


5 


0.1593 


0.837 


0.050 


6.20/5 


0.287 


4 


10 


6 


6 


0.1600 


0.863 


0.063 


5.80/5 


0.327 


4 


10 
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TABLE VI. Best ran^ 


fits to M(3Pi) 


for aniq 


= 0.025 using 


an axial vector current. 




i 


i 






error 


xVNdf 


conf. Ivl. 






1 


1 


0.1320 


1.695 


0.016 


15.91/2 


0.000 


9 


12 


2 


1 


0.1365 


1.543 


0.016 


13.63/2 


0.001 


9 


12 


2 


2 


0.1410 


1.430 


0.012 


1.47/4 


0.832 


7 


12 


3 


1 


0.1422 


1.337 


0.021 


21.97/2 


0.000 


9 


12 


3 


2 


0.1467 


1.223 


0.014 


1.38/4 


0.847 


7 


12 


3 


3 


0.1525 


1.000 


0.018 


1.82/4 


0.768 


7 


12 


4 


1 


0.1442 


1.317 


0.017 


20.52/4 


0.000 


7 


12 


4 


2 


0.1487 


1.145 


0.015 


2.06/4 


0.725 


7 


12 


4 


3 


0.1545 


0.911 


0.021 


2.13/4 


0.712 


7 


12 


4 


4 


0.1565 


0.816 


0.026 


2.01/4 


0.734 


7 


12 


5 


1 


0.1452 


1.275 


0.020 


12.19/4 


0.016 


7 


12 


5 


2 


0.1497 


1.099 


0.018 


2.59/4 


0.629 


7 


12 


5 


3 


0.1555 


0.861 


0.024 


2.50/4 


0.645 


7 


12 


5 


4 


0.1575 


0.762 


0.029 


2.39/4 


0.665 


7 


12 


5 


5 


0.1585 


0.706 


0.036 


3.04/4 


0.551 


7 


12 


6 


1 


0.1460 


1.236 


0.025 


5.82/4 


0.213 


7 


12 


6 


2 


0.1505 


1.061 


0.023 


4.26/4 


0.372 


7 


12 


6 


3 


0.1562 


0.822 


0.029 


3.21/4 


0.524 


7 


12 


6 


4 


0.1583 


0.720 


0.036 


3.14/4 


0.534 


7 


12 


6 


5 


0.1593 


0.772 


0.081 


0.54/3 


0.910 


8 


12 


6 


6 


0.1600 


0.797 


0.148 


0.13/3 


0.988 


8 


12 
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TABLE V 


. Best ran^ 


^e fits to M(iPi) 


for aniq = 


0.010 using 


the cubic group 


state ^Pi. 




i 


J ( 




aM(iPi) 


error 


xVNdf 


conf. Ivl. 


t^ 




1 


1 


0.1320 


1.713 


0.012 


1.81/4 


0.770 


5 
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TABLE IV 


. Best ran^ 


fits to M(3Po) 


for aniq = 


0.010 using 


the cubic group 


state Ai. 




i 


j (k, + Kj)/2 


aM(3Po) 


error 


xVNdf 


conf. Ivl. 


U 




1 


1 
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5 
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TABLE III. 


Best ran^ 


fits to M(3Pi) 


for aniq = 


0.010 using 


the cubic group 


state Tia- 




i 


i (« 


''i ~\~ '^»i)/2 




error 


xVNdf 


conf. Ivl. 


U 




1 


1 
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TABLE II. 


Best ran^ 


fits to M(3P2) 


for aniq = 


0.010 using 


the cubic group 


state Ea- 




i 


i (« 


''i ~\~ '^»i)/2 


aM(3P2) 


error 


xVNdf 


conf. Ivl. 


U 




1 


1 


0.1320 


1.706 


0.013 


1.14/4 


0.888 


5 


10 
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0.1365 


1.537 


0.015 


0.25/4 
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0.043 


14.88/5 
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for the aniq = 0.010 dynamical fermion and the quenched data sets. The other dynamical 
fermion data set yielded consistent results, but with much larger errors. See Tables XXI 
and XXII for results. 

VI. CONCLUSIONS 

This simulation had a number of inadequacies. First, our P-wave data are much noisier 
than S-wave spectroscopy. Lattice simulations with light fermions need better interpolating 
fields for P-wave states. It is much easier to explore large classes of trial wave functions and 
to find optimal ones when one has one or more heavy or static quarks (cf. Refs. [29] and 
[30]), and so heavy Wilson simulations are presently just not competitive with simulations 
with one static or nonrelativistic quark from the point of view of statistical uncertainties. 
Second, heavy Wilson fermions exhibit known lattice artefacts which should be corrected 
in future simulations through the use of improved actions. NRQCD remains the method of 
choice for heavy quark systems. However, the large scaling violations between relativistic 
fermion simulations and NRQCD need to be understood. Nevertheless, it is possible to make 
a determination of the strong coupling constant with an uncertainty comparable with other 
recent measurements, and which does not have the associated uncertainty induced by the 
absence of sea quarks. 

The final result for the strong coupling constant from this work is a^^~^^(Af^) = 
0.111(6), or agL=3)(M^) = 0.265(32). 
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interpolate linearly between these lattice masses to find the hopping parameter for which 
aM^ = a(3.f GeV). The error in the lattice spacing gives us a range for Kcharm- 

For the aniq = 0.010 dynamical fermion data set, we find Kcharm = 0.134 it 0.002, while 
for aniq = 0.025 we find Kcharm = 0.129 it 0.003, and for the quenched data set Kcharm = 
0.128 ± 0.003. 

We can compute the masses of the D mesons by looking at our heavy-light meson states, 
extrapolating the light quark mass to zero, i.e. Kught — ^ '«crit, and interpolating the heavy 
quark mass to the charm mass, Kheavy — ^ '^charm- For the nj = 2 data we divide our six 
fiavors of valence quarks into heavy and light as follows: Kheavy £ '«2, '^s} = {0.1320, 
0.1410, 0.1525} and Kught £ {^-4^ n-s^ n-e} = {0.1565, 0.1585, 0.1600}. The quenched data 
set has two heavy kappas and three light ones: Kheavy £ {'«o,'«i} = {0.1300, 0.1450} and 
'flight £ {'«2, '«3, '^4} = {0.1520, 0.1540, 0.1550}. Using the jackknife procedure mentioned 
above, for each Kheavy separately, we use the best fit ranges to the three mesons with different 
light quarks and extrapolate to Kught = '«crit- Having done this for the values of Kheavy we 

can interpolate to Kheavy = '^charm- 

We must shift the meson mass as with the J/V', but since the D mesons have only one 
heavy quark we replace Eqn. (13) by 

Mod = fi — niia + m2a, (40) 

with the identities (14) and (15). The results are summarized in Tables XVIII, XIX, and 
XX. Our results are consistently lower than those from experiment [21,28]. Using the 
cubic group representations we found complete degeneracy among the P-wave states, in 
disagreement with experiment. 

We also calculate the mass differences 

M{D*) - M{D''), (41) 

and 

M(i)i) - i[M(i)°) + 3M(i)*)], (42) 
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{Mr = 1.777 GeV) 



0.265(32), 



(37) 



which is in loose agreement with CLEO. 

Following the same procedure using the data with sea quarks of mass aniq = 0.025 rather 
than 0.010 gives the following numbers: 



Since the data seems to be less noisy for the aniq = 0.010 data than the aniq = 0.025 
data, we will use the former in quoting our final results, (36) and (37). 



Although there is already evidence that Wilson quarks do not calculate the hyperfine 
structure correctly and it is known that heavy Wilson fermions have an anomalously small 
lattice artifact magnetic moment [27], we attempt to calculate the masses of the Z = and 
L = 1 D-meson states. 

In Section IV A we estimated the hopping parameter of the charm quark, Kcharm, in 
order to calculate the S-P mass splitting of charmonium. We could only guess roughly at 
Kcharm siucc wc did uot kuow the lattice spacing. However, the S-P splitting is very weakly 
dependent on the quark mass, so we could determine it without knowing Kcharm- On the 
other hand, while the mass splitting is independent of quark mass, the masses of the charm 
states are not. Thus having found the lattice spacing, we can find Kcharm accurately and 
compute the mass spectrum of the charmed meson system. 

As we mentioned in Sec. IV A we must shift the lattice meson mass using Eqn. (13). In 
what follows we perform a jackknife analysis, dividing our 100 configurations into 10 sets of 
90 sequential lattices for the dynamical fermion simulations and dividing our 79 quenched 
lattices into 7 sets of 68. We find the best fits to the vector meson at our three lowest values 
of the hopping parameter and correct the heavy quark mass using Eqns. (13)-(15). Then we 
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V. CHARM SPECTROSCOPY 
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FIG. 10. Recent measurements of the strong coupling constant at the Z-mass. The square 
marks our result. Diamonds represent other lattice calculations, and crosses are experimental 
results. 

We now compare our number to other lattice determinations of the strong coupling 
constant and to experiment. It has become conventional to make these comparisons at the 
Z pole. Fig. 10 shows a compilation of recent lattice and experimental results, from Ref. 
[24]. 

It is interesting that, although the lattice spacings using Wilson quarks are significantly 
lower than those found on the same configurations using non-relativistic quarks, the deter- 
minations of the strong coupling are in good agreement [4] . 

The world average [25] for the strong coupling at the Z boson mass is ajjg{Mz) = 
0.117(5). Our number (36) is in agreement with this average. 

Running to the Z mass somewhat artificially compresses all the uncertainties of a \ow-Q^ 
calculation (the = fixed point is at infinite Q), and it might be more revealing to compare 
the lattice prediction to an experimental measurement at a low energy scale. The CFEO 
collaboration [26] gives the strong coupling at the r lepton mass as aj^^MT-) = 0.309(24). 



By running the MS coupling with nj = 3 down to this scale we find a lattice prediction of 
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modified minimal subtraction scheme, aj^. The two couplings are related perturbatively 
by 

«M5(9) = «y(9e'/')(l + ^ + 0K)). (29) 

Inserting (28) into the above expression gives 

ag^=^)(2.82 GeV) = 0.218 ± 0.021 ± 0.007, (30) 

where the first error is propagated from the error in ay, and the second is a systematic error 
of ay assumed from the perturbative expression (29). 



Once we have converted to the MS subtraction scheme, we can include the 3-loop term 
in the /3-function (22), (see [23]): 

- ;h C^) [lo, - 1] - - m(i - :^)|, (31) 



\ai{q), 

where 



A = 5 ^2857 - — „; + -ni) . and = (32) 

Thus our error from running the coupling with this /3-function should be of order a'*, which 



is always smaller than the present errors. 



Rodrigo and Santamaria [23] show that the heavy quark thresholds in MS prescription 
are at the mass of the quark. When we cross the threshold, we increment by one. Starting 
with 

ag^=^)(2.82 GeV) = 0.218 ± 0.022, (33) 

we run down to the charm quark mass with three fiavors, and change to four. Then we can 
run to the bottom quark mass and change from four fiavors to five. Finally, we can run the 
coupling up to the Z mass with five fiavors. The results follows: 

a^^^'^\q = 1.5 GeV ~ m,) = 0.287(38) (34) 
a^^^'^\q = 4.5 GeV ~ m^) = 0.195(17) (35) 



a^^^\q = 91.2 GeV = Mz) = 0.111(6). (36) 
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we must run one of the couplings to the scale at which the other is known using the two-loop 
/3-function (22). 

We need only one value for ci'y^ , so for the time being we shall use the lattice spacing 
from the am^ = 0.010 staggered fermion data. Thus we have the strong coupling measured 
with either zero or two flavors of sea quarks: 

a[r'^=°^(6.14 GeV) = 0.152 ± 0.005 ± 0.005 

a[r'^=^^(6.48 GeV) = 0.179 ± 0.004 ± 0.010, (26) 

where the first errors quoted are propagated from the uncertainty in the lattice spacing, 
and the second errors are our estimation of the systematic uncertainties resulting from the 
perturbative calculations. We choose to run the fiavor coupling to 6.48 GeV, giving 

0.152 ± 0.006 for = 

ay(6.48 GeV) = <^ (27) 

I 0.179 ± 0.011 for nj = 2 

Here the statistical and systematic uncertainties have been combined in quadrature. If we 
had chosen to perform the extrapolation to three fiavors at 6.14 GeV instead, there would 
be no change in our determination of ay at any scale. 

We may linearly extrapolate either 1/a or a from zero to three fiavors. (In first order 
1 /a depends linearly on n J.) Extrapolating in 1/a gives cty ''^ (6.48 GeV) = 0.198ib0.020. 
However, linear extrapolation in a has smaller errors, 

alr''=^V-48 GeV) = 0.194 ± 0.017, (28) 

and we will use it below. 

D. Lattice Prescription to Minimal Subtraction Schemes 

In order to compare our result with experimental determinations of the strong coupling 
constant, we convert from our lattice definition of ay (28) to the coupling defined in the 
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0{ay) is a source of systematic uncertainty. Thus we invert (9) neglecting higher order 
terms and estimate our systematic error to be (1.19 + 0.07nj)^ay. 

Using the procedure outhned above, we find that, for our dynamical fermion simulations 
where (Tr [/piaq)/3 = 0.5650 and 0.5644 for aniq = 0.010 and 0.025 respectively, 

. 0.179 ± 0.010 for am, = 0.010 

ai;^-^)(3.41/a) = ' (21) 

I 0.179 ± 0.010 for am, = 0.025 

where the errors quoted are systematic. We do not quote an uncertainty for the trace of the 
plaquette because it is much smaller than our other errors. 

Since we are using the mass splitting in the charmonium sector, three fiavors of sea 
quarks contribute to the theory. There are two methods we know of to convert ay from 
nj = 2 to nj = 3. In the first we run ay down to some low momentum scale (~ 500 — 1000 
MeV) with two fiavors using the two-loop /3-function, 

aiq) = aiiq) + aiiiq) + Oia^), (22) 



where 



1 + a{qo)Pot 

au{q) = -[a,{q)]\log(^), (24) 



wi 



ith the definitions 



t = ^log(pj, /3o = ll-^nj, /3i = 102-^nj, and W = (25) 
47r \qoy 3 3 47rpo 

However, this method is inadequate because we do not know the scale of the strange quark 
threshold, and because our perturbative expression (22) fails as ay becomes large. 

The second method is to use our quenched data to give us the strong coupling for nj = 0. 
From (Tr [/piaq)/3 = 0.59367 we find that Jv^^^^SAl/a) = 0.152(4). Since the inverse 
coupling is nearly linear in nj, we can extrapolate to nj = 3. However, since the lattice 
spacings are different for the dynamical fermion simulations and the quenched simulation, 
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degenerate using Wilson quarks. However this is not the case in the physical world; i.e. 
using the ^Pi state to fix the scale rather than the ^Pi state changes our determination of 
the lattice spacing. Thus we must assign a value to the systematic uncertainty to of 
about 100 MeV which we add in quadrature to a statistical error of about 50 MeV to find 

a-^ = 1900 ± 110 MeV for arrig = 0.010. (18) 

Similarly, we quote 

a-^ = 1660 ± 150 MeV for arrig = 0.025. (19) 

One puzzling feature of the calculation is the very different lattice spacings found here 
using heavy Wilson quarks and the lattice spacing of the NRQCD group. The numbers are 
shown in Table XVI, along with all other lattice spacings extracted from this data set of 
which we are aware: in the table S and W label staggered and Wilson valence quarks, the 
"force" is from the string tension [22], and the zero quark mass line is from an extrapolation 
when it is available. 

Presently we have calculated only the ^Pi — S mass difference on our quenched configu- 
rations (see Table XVII). The lattice spacing given by that measurement is 

a-^ = 1800 ± 180 MeV for arrig = oo. (20) 

We should point out that for this mass difference, it was necessary to fit the axial-vector 
propagator to a different range than the pseudoscalar and vector propagators. The latter 
two take a longer time to "saturate" the signal, and the axial-vector propagator becomes 
noise dominated after some time. The ranges used are those indicated in Tables XI-XIII. 

C. The Lattice Calculation of 

The expectation value of the plaquette determines the coupling on the lattice at mo- 
mentum scale 3.41/a. However, the fact that we do not know the coefficient in front of the 
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propagators together in order to compensate for correlations due to using the same config- 
urations for all three quantities. The results for the quarkonium [nq = Hq) pairs are listed 
in Table XIV for dynamical staggered fermion mass = 0.010, Table XV for fermion mass = 
0.025, and in Table XVII for the quenched approximation. 

Since we cannot know a priori the value of Kcharm, and since we assume AM to be 
insensitive to n for heavy mesons, we average AM for the two heaviest quarkonium pairs 
and set them equal to the physical value Mx - (M^ + 3M^J/4 = 442 MeV [21] to find the 
lattice spacing. Again, one should refer to Tables XIV, XV, and XVII for results. 

Instead of using an interpolating field to couple to the ^Pi state, we can couple to the 
full P-wave system by using the corresponding lattice multiplet. Thus, we can also compute 
the mass differences 

^M{^Pr-S) = MCP^)-^[MCSo) + mCSi)], (16) 

and the true "spin-averaged S-P mass splitting", 
AM(P-:^) = ^[5MfP2) + 3MfPi) + MfPo) + 3M(Vi)] - i[M05'o) + 3MfS'i)] 

(17) 

To do this, we take the best operator corresponding to each J^*^ value and in the case of 
(16) perform a three-propagator correlated fit, or in the case of (17) a jackknife fit. So far 
this has been done for the dynamical fermion simulations only. The lattice mass differences 
are presented in Tables XIV and XV. 

B. Lattice Spacings 

Tables XIV and XV give a summary of the inverse lattice spacings obtained from the 
S-P mass splittings. Each mass splitting gives a lattice spacing, thus we must decide which 
value to use to set the scale. For the aniq = 0.010 data set, we note two things: the lattice 
values for the P — S difference are the same which implies that the ^Pi and ^Pi states are 
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Since we are going to fix our measured lattice splitting to the mass splitting in charmo- 
nium, we need to have a rough idea of Kcharm- In previous work on these configurations [8], 
the lattice spacing was determined by fixing the rho to its physical value, yielding 1/a = 2140 
and 2000 MeV for sea quark masses 0.010 and 0.025 respectively. To find the approximate 
value for Kcharm we use these lattice spacings to find the approximate k, for which the J/t/j 
mass is at its physical value, i.e. M^a ~ 3.1/2.0 = 1.55. Also, since the charm quark mass 
is of the same order as the inverse lattice spacing [nic <^ «~^), the dispersion relation for 
Wilson fermions is distorted: 

E{k) = m, + ^ + ... (12) 
2m2 

with nil ^ 1712- To account for this we must shift our lattice mass [19,20], using 

Af^a = fi — 2mia + 2m2a, (13) 

where, with k = K/(8/€crit), 

/\ — Qk \ 

niia = log y 2k J' ^^^^ 

and 

exp(mia) sinh(mia) 

= — — . (15) 

1 + smh(mia) 

'«crit is estimated from [8] to be 0.161 for both values of sea quarks, and 0.157 for the quenched 
data. This gives us the estimate that for the dynamical fermion data set Kcharm is between 
0.1320 and 0.1410, the values of our two heaviest fiavors of valence Wilson quarks, and just 
below 0.1300 for the quenched data set. Of course, this was done with a rough guess for the 
lattice spacing, which is indeed the quantity we are attempting to calculate. However, since 
we argue that the S-P splittings are insensitive to the quark mass in the heavy systems, we 
need only a rough idea of the value of Kcharm- 

Therefore, for a given pair of k, values, we fit the correlators using pseudoscalar, vector, 
and axial vector currents to (11). When calculating the mass difference we fit the three 
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Once the lattice spacing and the coupling constant are fixed one may convert the coupling 
constant to other schemes and run it to other desired scales. 

Of course, this calculation is incomplete. It depends on the empirical observation that 
the coupling constant defined through the potential provides a good perturbative expan- 
sion parameter from fairly large lattice spacing on down. A completely nonperturbative 
alternative is the calculation of the strong coupling constant by Liischer, Weisz, Wolff, and 
collaborators [16]. Here the coupling constant is defined via the response of a lattice system 
to an external color electric field, and a wide range of physical scales can be covered by a 
series of steps each of which involves only a small change of physical scale. From a phe- 
nomenological point of view the two methods produce equivalent results for ajjg{Mz) from 
quenched simulations. This program has not been carried out for full QCD yet [17]. 



To set the scale we use the mass difference between P-wave states [L = 1) and S-wave 
states [L = 0) in the charmonium sector. We choose this physical quantity, rather than 
the mass of the rho meson for example, for two reasons. First, the S-P splitting is fairly 
insensitive to the quark masses. Second, in the laboratory the widths of the charmonium 
states are much narrower than those of the light mesons. Therefore, it is realistic to think 
that the S-P mass splitting is less sensitive to lattice artifacts than are light hadron masses. 

Using currents to couple to the pseudoscalar, vector, and axial vector mesons, we can 
compute the mass difference 



We perform a six parameter correlated fit to the three propagators in order to calculate this 
mass difference [18] using a fit function for each propagator of 



A. Fixing the Lattice Spacing 
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FIG. 9. "^Pi — ^Pi meson mass difference from the arrig = 0.01 data set as a function of the 
average quark hopping parameter at k = (Kg + Kg)/2. 



IV. THE STRONG COUPLING CONSTANT 

The lattice allows a determination of the strong coupling constant through the mea- 
surement of two physical and nonperturbative quantities. The mass splitting between two 
states, computed on the lattice, gives a length scale a, the lattice spacing. The plaquette (or 
any other short distance quantity) gives a coupling constant at a momentum scale inversely 
proportional to the lattice spacing. 

The method uses a definition of the strong coupling constant in terms of a physical 
observable [15]: that is, we define av{q) through the potential 

V{q) ^ .^ilvrovM. (8) 

Here C/ = 4/3 is a group theory factor, and q is the gluon momentum. With this definition, 
the perturbative expression for the logarithm of the trace of the plaquette is [12,15]: 



log ( ^Tr Upiag'j = yay(3.41/a) 1 - (1.19 + 0.07nj)ay + 0{ai 
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FIG. 7. ■^P2 — "^Pi meson mass difference from the arrig = 0.01 data set as a function of tlie 
average quark flopping parameter at k = (Kg + Kg)/2. 
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FIG. 8. "^Pi — "^Po meson mass difference from tlie arrig = 0.01 data set as a function of tlie 
average quark flopping parameter at k = (Kg + Kg)/2. 
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FIG. 6. Mass difference between two different representations of the cubic group corresponding 
to the meson, from the anig = 0.010 data set as a function of the average quark hopping 
parameter at k = (Kg + Kg)/2. 

All lattice data are correlated, and so to look for fine structure by simply subtracting 
the two masses overestimates uncertainties. Instead, we perform a correlated fit of the two 
propagators to two masses and extract an uncertainty from the error matrix. We present 
pictures of these differences in Figs. 7-9 (for the — "^Pi, ^Pi — "^Pq, and ^Pi — ^Pi mass 
differences, respectively). For the heaviest quark masses we are unable to see any fine 
structure splitting within our statistical uncertainty. 
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ignore them in the rest of the paper. 

Again, we compare effective mass fits and fits to a range of points. In this case we have 
several operators for each combination of spin and parity. Generally one or two operators 
provide a superior signal compared to the other ones. 

We will restrict our presentation of tabular masses to the following states: for the 

state we use one of the two E states, which we will refer to as Ea'. + O ). Here 

the first subscript refers to the "spin" content with it standing for ^f± and for 73, while 
the second subscript refers to the orbital angular momentum content with it standing for 
xrLiy and for 2;, up to normalization factors. For "^Pi we use the Tia state, ;^((9o+ — O+o), 

for ^Pq the Ai state, -^{Oj^ + O |_ + Oqo), and for the ^Pi we use the cubic group ^Pi 

state, 05+, where the 5 refers, of course, to 75. Note, again, that all the states used couple 
to the orbital x -\- iy source. The lattice masses for best fits for each different spin-parity 
combination are listed in Tables II-V for the aniq = 0.010 data and in Tables VII-X for the 
aniq = 0.025 data. 

We found that the different representations of the cubic group corresponding to the same 
angular momentum states are degenerate: the lattice does not break rotational symmetry so 
badly that we can observe it in spectroscopy. Fig. 6 illustrates this point: it is a correlated 
fit to the mass difference of the Ea and (:^(^o+ + O+0)) states, both of which correspond 
to ^P2 continuum states. 
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B. Cubic Group Representations 



On the lattice, instead of the angular momentum group, the cubic group is the symmetry 
group according to which states form irreducible representations. The "current" operators 
previously discussed have "orbital cubic group" representation Ai, corresponding to "orbital 
angular momentum" representation L = 0. The "spin cubic group" representations are Ai 
for the pseudoscalar and Ti for vector and axial vector currents. 

We can also build the equivalent of "orbital angular momentum" Z = 1, i.e., P-wave, 
states on the lattice. They are in the cubic group representation Ti and have wave functions 
xf[r)^ yf{r) and zf[r) with /(r) a function depending only on the radius. In our case 
/(r) is taken as a Gaussian. Equivalently, the wave functions can be taken as [x + iy)f[r)^ 
[x — iy)f{r) and zf[r). On a finite lattice, in order to satisfy periodic boundary conditions 
in the spatial directions, we use sin(27ra;/;/iv/;) instead of Xk in the wave functions. 

This "orbital cubic group" representation can be combined with the two different "spin 
cubic group" representations to get 

ri(x)Ai = ri, (6) 

Ti (X) Ti = Ai + Ti + ^ + r2 (7) 

The first line corresponds to a continuum ^Pi state and the second to continuum ^Pq, ^Pi 
and ^P2 states. The lattice E and T2 states combine in the continuum limit to become the 
^P2 states. 

In the numerical measurements we use one "orbital P-wave" (cubic group Ti) source, 
the X -\- iy one. We use all spin combinations at the source and sink and all "orbital cubic 
group Ti" wave functions at the sink. With this we can build 7 of the 9 possible correlation 
functions of the ^P states. One each of the three Ti and T2 states are inaccessible with our 
single X -\- iy source. In addition one can insert extra 70 's at either sink and/or source. We 
implemented only the variant with 70 's at both source and sink. However, the result turned 
out to be somewhat noisier than without these 70 factors in the operators and thus we will 
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and select the best one, in our case the one which maximizes the confidence level times the 
number of degrees of freedom divided by the statistical error of the mass. (This selection 
criterion has been used in our previous studies [8,12].) 

Effective mass plots and best fits to the dynamical staggered fermion data for the S- 
wave states have been published in a previous work [8]. In Figures 1 and 2 we show effective 
masses for the ^Pi state for aniq = 0.010 and 0.025. Note that, unlike the S-wave states, we 
can only fit out to distances t ~ 10 since we lose the signal in noise. The best range fits are 
listed in Tables I-V and VI-X. 

Results from the quenched simulations are presented here as follows: effective mass plots 
for the pseudoscalar, vector, and axial-vector currents are shown in Figures 3-5, and our 
best range fits are presented in Tables XI-XIII. 
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FIG. 1. Effective masses for the "^Pi state using Ai = for aniq = 0.010 at 
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which for large t reduces to a single decaying exponential 



C(t) ~ (O|0|l) 



where m is the mass of the lightest particle that couples to O. We employ two types of 
operators to couple to the P-wave states: operators whose spatial dependence is symmetric 
and which couple to the desired state through the quarks' spins ("current" operators), and 
operators which form specific representations of the cubic group. 



A. "Current" Operators 

The first class of operators we use are the analogs of currents. In this work we use the 
pseudoscalar, vector, and axial- vector currents given respectively by 

P = V-TsV-, (3) 
K- = V'7.V', (4) 
A = V'7j75V'- (5) 

These currents have respective quantum numbers J^*^ as follows: '", 1 , and 1"'""'". In 
spectroscopic notation, ^'^~^^Lj, which we shall use hereafter, these states are the ^5*0, "^Si, 
and "^Pi, respectively. 

For the present work we consider correlators which have "shell" sources and "point" sinks. 
That is, we measure the source of the correlator over a spatial volume with a Gaussian 
distribution centered at a point, and the sink at a single point. We also have data with 
"shell" sinks, from which the masses are consistent with, but noisier than, the point sink 
data. In either case the sink state is projected onto k = 0. 

One can examine the quality of the signal by calculating "effective" masses: fits to the 
correlator over just two neighboring points. If the assumption that the current is coupling 
to just one particle state is correct, then we should see a plateau in the effective mass. Once 
we establish an approximate range over which to fit, we fit every possible range therein 
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HEMCGC collaboration with two flavors of dynamical staggered quarks [8]. The configura- 
tions were generated using the hybrid molecular dynamics (HMD) algorithm [9]. The size of 
the lattices is 16"^ X 32, the lattice coupling is /3 = 5.6, and the dynamical quark masses are 
aniq = 0.010 and 0.025. Periodic boundary conditions were used in all four directions of the 
lattice. The total simulation length was 2000 simulation time units (with the normalization 
of ref. [10]) at each quark mass value. We analyzed lattices spaced by 20 HMD time units, 
for a total of 100 lattices at each mass value. 

The spectroscopy was computed with six values of the Wilson quark hopping parameter: 
K = 0.1600, 0.1585, 0.1565, 0.1525, 0.1410, and 0.1320. The first three values are rather light 
quarks (the pseudoscalar mass in lattice units ranges from about 0.25 to 0.45), and the other 
three values correspond to heavy quarks (pseudoscalar mass from 0.65 to 1.5). Our inversion 
technique is conjugate gradient with preconditioning via incomplete lower-upper (ILU) de- 
composition by checkerboards [11]. For more details about the dynamical staggered fermion 
simulations see ref. [8,12]. Since we use sources for the propagators which are extended in 
space, we fix gauge to lattice Coulomb gauge using an overrelaxation method [13]. 

The quenched configurations were generated using an updating algorithm which treats 
each link to a combination of four microcanonical overrelaxed hits on each of three SU(2) 
subgroups, followed by one Kennedy-Pendleton [14] quasi-heat bath update. The lattices are 
16"^ X 48 at a lattice coupling /3 = 6.0. Lattices were recorded for analysis every 200 sweeps, 
and we acquired a total of 79 lattices. Here, the spectroscopy was computed with five values 
of the Wilson hopping parameter: k = 0.1550, 0.1540, 0.1520, 0.1450, and 0.1300. Our 
boundary conditions, gauge fixing, and inversion technique are identical to those mentioned 



above. 



III. P-WAVE SPECTROSCOPY 



We extract masses from correlation functions of some operator (9, 



at) 



(0|O(t)O(0)|0). 
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The Fermilab group [3] was the first to use the S-P mass sphtting of heavy mesons to set 
the scale, which allows one to run the coupling calculated from the plaquette to any scale. 
However, since their calculation was in the quenched approximation, they had to correct for 
the presence of dynamical fermions. Our simulations are done in the presence of dynamical 
staggered fermions, so we avoid this extrapolation. Similar calculations have been done 
in two-fiavor QCD [7] and using non-relativistic quarks [4]. Further details regarding the 
simulations are given in Section IL 

We use interpolating fields which are generalizations of local currents to couple to the 
two S-wave states and to one of the P-wave states (the ^Pi). However, in order to create all 
of the P-wave states, we employ several representations of the lattice cubic group. Both of 
these procedures are described in Section HL 

In Section IV we describe in detail our methods for extracting the strong coupling from 
measurements of the plaquette, for using zero fiavor and two fiavor calculations of to 
extrapolate to three fiavors, for changing from our lattice definition of to the modified 
minimal subtraction definition, and for perturbatively running the coupling to any mo- 
mentum scale. Particular attention is given to possible sources of systematic errors in our 
calculation. 

Finally, in Section V we present lattice calculations of the masses of charmed mesons. 
Our computations confirm that the Wilson action does not generate the hyperfine structure 
correctly. 

II. THE SIMULATIONS 

Our dynamical fermion simulations were carried out on the Connection Machine CM-2 
at the Supercomputing Computations Research Institute at Florida State University, and 
our quenched simulations on the Paragons at the San Diego Supercomputer Center and at 
Indiana University. 

For the dynamical simulation we used the ensemble of configurations generated by the 



3 



I. INTRODUCTION 



This paper presents results of P-wave meson masses using Wilson quarks on both 
quenched gauge configurations and configurations which include light dynamical quarks. 
Here we focus on systems containing one or two heavy quarks. 

Lattice calculations of P-wave systems go back only a few years. Simulations with stag- 
gered fermions regularly measure the masses of states which are odd-parity partners of the 
ordinary ground state mesons. The earliest calculations with Wilson fermions of which we 
are aware were done by the APE collaboration [1] and then in 1992 one of us [2] presented 
incomplete calculations of the whole P-wave multiplet, which hinted at the existence of fine 
structure splitting in charmonium in qualitative agreement with experiment. Also in 1992 
the Fermilab group performed P-wave spectroscopy for heavy quark systems with Wilson 
fermions and with improved Wilson fermions [3]. Since then very precise calculations in 
heavy quark systems have been done using simulations with nonrelativistic quarks by [4-6], 
with much higher accuracy than Wilson fermions permit, for the same amount of computer 
time. Indeed, the calculations we present do not show statistically significant fine structure 
splitting at our heaviest quark masses. 

Why are more calculations with Wilson fermions being presented? There are several 
reasons. First of all, the spectroscopy of states with light quarks demands a lattice imple- 
mentation of relativistic quarks, so in D meson spectroscopy (for example) at least one of 
the quarks must be relativistic. The heavy quark may be treated on the lattice as static, 
or nonrelativistic, or relativistic but heavy. Fach prescription has its own set of systematic 
errors, and so it is important to do the calculation in all three ways on the lattice to under- 
stand them. The calculations done in this paper for charm spectroscopy complement ones 
using the same gauge configurations in a companion work [6] using nonrelativistic heavy 
quarks. 

Second, even without particularly accurate P-wave mass measurements we can address 
a topical physics problem: the extraction of the strong coupling constant from lattice QCD. 
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Abstract 

In this work we present lattice calculations of the masses of P-wave mesons 
using Monte Carlo simulations. Our valence fermions are defined by the 
Wilson action. Our gauge fields are generated with both dynamical staggered 
fermions at a lattice coupling (3 = = 5.6 for sea quark masses of anig = 
0.010 and 0.025, and in the quenched approximation at /3 = 6.0. We present 
results for charm and charmonium spectroscopy and use them to compute the 
strong coupling constant a^. We compare our results to those of other recent 
lattice calculations and experiments. 
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